The Intrinsic Spin Hall Conductivity in a Generalized Rashba Model 
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We calculate the intrinsic spin Hall conductivity <r sH of a two-dimensional electron system within 
a generalized Rashba model, showing that it is, in general, finite and model-dependent. Considering 
arbitrary band dispersion, we find that cr sH in the presence of the linear-in-momentum spin-orbit 
coupling of the Rashba form does not vanish in the presence of impurities except for the precisely 
parabolic spectrum. We show, using the linear response Kubo formalism, how the exact cancellation 
happens for the quadratic dispersion, and why it does not occur in general. We derive a simple 
quasiclassical formula for <r sH in terms of the Fermi momenta for the two electron chiralities, and 
find that <j sH is in general of the order of the squared strength of the Rashba term. 
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tie spin Man ettect Q 
ject 0, E| or a rather old one Q depending on one's 
perspective. Much recent interest has focused on the 
intrinsic spin Hall effect (ISHE) where a non-magnetic 
system (e.g. a III-V semiconductor) spontaneously ex- 
hibits a bulk spin Hall conductivity er sIi (i.e., a spin 
conductivity transverse to the direction of an external 
electric field in the absence of any applied magnetic 
field) arising entirely from the spin-orbit coupling ef- 
fects in the underlying one-particle band structure. The 
ISHE is to be distinguished from its better-known coun- 
terpart, the extrinsic SHE (ESHE) @, 0, predicted 
by Dyakonov and Perel Q], which arises due to spin- 
correlated asymmetry in impurity scattering. Recent 
claims of experimental observations of both ESHE 0, 0] 
and ISHE (in semiconductor structures) make the 
subject particularly interesting since there is controversy 
about the existence of ISHE in various theoretical models 

mEEQHHHIIllIllIllIIl] In particular, 
within the extensively studied Rashba model [22l |. aris- 
ing from the spatial inversion asymmetry in 2D systems, 
the intrinsic cr sIi has been claimed to be a large univer- 
sal constant 3] e/8iT, where e is the electron charge, but 
subsequently it was found that an arbitrarily small con- 
centration of impurities in the Born approximation (i.e., 
when one neglects all spin-related asymmetry in scat- 
tering leading to ESHE) introduces vertex renormaliza- 
tion exactly canceling 0, O 0, 01 the universal 
ISHE predicted in the Rashba model in [3j. (The same 
cancellation also occurs for the 2D Dresselhaus model 
for the linear-in-momentum spin-orbit coupling.) For 
other models of spin-orbit coupling, e.g. the 3D Dres- 
selhaus model 12(1 o r the Luttinger model for the va- 
lence band holes [2l(, the ISHE has been theoretically 
calculated to be finite and model-dependent. Since the 
Rashba model is by far the most widely studied model 
for the spin-orbit coupling in 2D semiconductor struc- 
tures and since impurities are crucial for the validity of 
the linear response theory, the precise vanishing of the 
ISHE in the Rashba model has led to substantial con- 
fusion about whether the ISHE, as opposed to ESHE 



010001 which is alwa ys p resent, is ever finite in the 
Rashba model H2lli3T . il 4 llfil | and perhaps even in all 
models [lflj . 

We show in this Letter that the precise vanishing of 
the 2D ISHE in the Rashba model is limited to quasi- 
classical calculations for the usually assumed quadratic 
unperturbed band spectrum 



eo(p) =p 2 /2m. 



(1) 



In fact, the cancellation is accidental, and slight changes 
in the model (i.e., using a band dispersion different from 
the usual parabolic band dispersion) give rise to a fi- 
nite ISHE even in the Rashba model. While the vertex 
corrections do decrease the magnitude of intrinsic er sH , 
they contain an integral over momenta that depends on 
band dispersion eo(p), and cr sIi is finite for a general (i.e. 
non-parabolic) band dispersion, albeit explicitly model- 
dependent. Our finding of a finite 2D ISHE in the gener- 
alized Rashba model is of particular significance in estab- 
lishing a matter of principle, refuting the earlier findings 
of either the large universal ISHE or the vanishing ISHE. 

In the presence of spin-orbit coupling in a Hamiltonian, 
a spin current of the from 
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should appear in the kinetic equation for the spin density. 
This quantity leads to spin accumulation at the boundary 
of a sample, limited by spin diffusion, when an electric 
current is passed in the absence of magnetic field (the 
spin Hall effect). The appearance of such a current was 
predicted for a band structure of III-V semiconductors 
Q and for general 2D Rashba coupling Q , and has been 
a subject of a great deal of recent activity pj. 

We derive in this Letter the spin Hall conductivity 
for a general spectrum, which has a simple quasiclas- 
sical limit. First we demonstrate how the cancellation 
for the parabolic disp ersion, already noted in the liter- 
ature (ill H2I Il3l Q, ^^|, takes place and also why it 
is not universal but explicitly depends on the spectrum. 
We then provide the general quasiclassical limit for the 
ISHE conductivity. 
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Consider the two-dimensional electron Hamiltonian 



H = eo(p) + az(<7 x p) 



(3) 



comprising an unperturbed part eo(p) and the Rashba 
coupling (parameterized by the strength a) of spin and 
orbital movement in the plane defined by the normal z. 
From the velocity operator is 



Vi(p) = d Pt H = e' (p)pi + ae zji <jj 
The spin-velocity is usually defined as 



'(P) = 



(4) 



(5) 



In the framework of the linear response theory the spin 
Hall conductivity can be calculated as a bubble (Fig. 
with the vector vertices j = ev and j z = \w z . 

To emphasize the cancellation we break up the spin- 
Hall conductivity in two parts (the two terms in paren- 
theses) 



x(t$(p) + iy(w,n))G(p,w + fi)]. (6) 

The first term is the bare bubble in Fig. Q](a), the sec- 
ond term is the result of the summation of the ladder 
series in Fig. (b). The terms of the ladder series do 
not vanish after averaging over p because the Green's 
function G(p, lu) for the Hamiltonian itself depends 
on the direction of momentum. 

The contribution from the first (bare) term in JfjjJ is 
erfj 1 = e x ijOQ, where the quantity crg H is given by an 
integral that depends on the spectrum £o(p): 



where 



(7) 



(«) 



Here G u , rf (p,w) = (iuj - £ M ,d(p))~\ where £ Uid = £o(p)T 
ap — ep and S = u> + sgnw, is the Green's functions 
for the chiral eigenstates of the Hamiltonian . For 
brevity we omit the arguments of the Green's function, 
and denote the ones with w + Q by the subscript +. 

Taking the integral over momenta in in the £- 
approximation with quadratic spectrum we find (see 
below for details) 



(9) 



1 + 4A 2 r 2 ' 



where 2A = 2apF is the value of the spin-orbital energy 
gap for the two chiralities at the Fermi level. In the clean 
limit t — ► oo, Eq. gives to the universal constant 



e/87r Q independent of the strength of the spin-orbit 
interaction a. The physical meaning of this constant can 
be easily seen if one carries out the frequency summation 
in CTo H first. Then in the clean limit one gets 



(10) 



which at T = singles out the annulus pd < p < p u , pop- 
ulated by the electrons with the chirality u only, where 



p U: d = \l m 2 a 2 + 2mep ± ma are the Fermi momenta for 
the two chiral eigenstates for the spectrum JQ: 



7o H (0) = 



Pu pdp £q(p) e Pu 



I'd 



2ir (2ap) 2 8tt 2ma 



8tt' 



(11) 



In (0 the ladder summation is carried out in the form 
of the renormalization of the spin-velocity vertex v z — > 
v z + V z . The result for the summation of the ladder 
corrections is of course independent of which of the two 
vertices, the current ev or the spin-current ^v z , in 
is renormalized. Straightforward calculations (shown at 
the end of the paper) yield the ladder correction e Z ycr| , 
where 



aB[e ] - gjep] 
2-kvqt - B[e ] ' 



(12) 



Here 



B[eo] 
C[eo] 



lE P (Gu + + G d+ )(G u + G d ) (13) 
lJ2 P ^(p)[Gu + G u -G d+ G d }. (14) 



Except for the averaging over direction p of momentum, 
we did not perform momentum integration in the quan- 
tities -A[eo], B[eo] and C[eo]- Summing up (Tg 11 and crf^, 
we see that 



sH 



-tVa[ £0j 



olI-kvqt — C[eo] 



(15) 



An exact cancellation takes place only when C[eo] = 
a2-n:vQT. This is the central result of the paper. Such 
a cancellation is indeed what happens for the quadratic 
spectrum (0. 

In general, keeping eo(p) for band dispersion in the 
Green's functions, we get in the ^-approximation 



C[eo] = 2ttuqt 



4^ n 



(16) 



Here v Ui d is the density of states on the Fermi surface for 
the chiralities u and d, vn = {v u + Vd)j2 by definition, 
and I — g(sgn(o; + ft) — sgna>) is a factor accounting 
for the requirement that the two poles in the Green's 
functions in the bubble lie in different half-planes (i.e., 
one Green's function is retarded and one advanced). If we 
suppose fl > 0, then this factor equals unity in the band 
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a) 




contributions: 



b) 




c) 



< = 




Figure 1: Bubble diagrams for the spin Hall conductivity: 
bare (a) and with vertex renormalized by the ladder series of 
impurity scattering (b). Equation on the latter (c). 



— fi < uj < and zero otherwise. Such a factor generally 
appears in the frequency summation in cr sIi , which means 
that all quantities under the frequency summation should 
be taken in the specified frequency band, where we may 
simply deem 7 = 1. For fl < the reasoning is similar. 

We note that the density per chirality in 2D irrespec- 
tive of the spectrum is n u = p^/4ir, hence dn u /dp u = 
p u /2-K. On the other hand, using that, by definition, e F = 
£o(j>u)-aPu, we get dn u /dp u = v u de F /dp u =u u (^ (p u )- 
a). For the d chirality, we get analogous expressions with 
the opposite sign before a. We finally get 



4z/n 



2 M 



Pu -Pd 



(17) 



We therefore reduced the expression that explicitly de- 
pended on the electron spectrum by a simpler one that 
still depends on the spectrum implicitly through the 
Fermi momenta p UtC t for the two chiralities. In this new 
form, Eq. II17J1 . however, we immediately see whether 
the exact cancellation for ISHE takes place for a given 
spectrum. For the quadratic dispersion Q p u — Pd — 
2ma — iani> , where Uq — m/2n, so that (flTI) gives a, 
and C[eo] = c&ttvqt exactly cancels out the spin Hall 
conductivity 11511 . 

Above we used the ^-approximation in its usual form, 
i.e. the prefactors before the Green's functions were 
taken at the Fermi level. It may be shown that for the 
quadratic spectrum the cancellation is exact in all 
orders of the quasiclassical approximation, i.e. if one ex- 
pands prefactors in powers of £, all the terms vanish. 

We can calculate the total spin Hall conductivity for an 
arbitrary spectrum in the quasiclassical ^-approximation. 
Using the expressions for B and ^4[eo] obtained in the end 
of the paper, we get for the sum of the bare and ladder 



Pu 



Pd 



8n \ Aitavo 



1 



(18) 



For a spectrum different from JQ, the cancellation in 
the expression above does not take place, e.g., if eo(p) = 
vop, then p U: d = e F / (vq ± a) and 



v = 



e F v 2 + a 2 



2ttv 2 



or 



(19) 



So that p u — pd = 47r^oa(«o — ct 2 )/{vQ + a 2 ), whence 
CT sH = -(e/^)2a 2 /{vl + a 2 ). 

It may be shown that the net <r sH is always at most 
of the order a 2 . Expanding p Ut d in ijlSjl in powers of 
a around the Fermi momentum p F of the general band 
dispersion £o(pf) = £f, we find that the lowest order 
expression always cancels the unity in the parentheses of 
Eq. II18J1 . and leaves the next order term, 0(a 2 ), as the 
leading ISHE contribution for small a. 

Below we provide details of the derivation of our results 
above. The Hamiltonian 10 can be diagonalized: 



U + 7iU = €q(p) — apa z 



by rotation in spin space: 



1 



U=±( 1 



(20) 



(21) 



where ipp is the angle between p (the direction of elec- 
tron momentum) and x. The eigenstates of the rotated 
Hamiltonian (l2f))) are spin up and spin down. In the 
original basis they correspond to the two chiralities (spin 
along and opposite to z x p). Columns of U give the 
eigenfunctions in the original basis. 

In the basis of the chiral eigenstates the Green's func- 
tion, which is a matrix in spin space, is diagonal with the 
elements G Ut d- In the original frame 



G(p,w) = U 



IP 



G u 
G d 

= \{G U + G d ) - \{G U - G d )e zl] a lPr (22) 



The renormalized vertex V 2 is found from an inhomo- 
geneous linear equation 



1 



27r^oT 



G(P, «) (v z (p) + V 2 ) G(p, w + O), 



(23) 

represented by the diagrammatic equation in Fig. (c) . 
Here 1/27tz/ot is the impurity line in the Born approxi- 
mation. The inhomogeneous term, being a spin matrix, 
has only off-diagonal elements (i = cc, y) 

J2 G(P, (p)G(p, w + fi) = o t A[ea] . (24) 
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Hence V 2 only has the same components (that depend 
on oj, fl and not on p): 



V 2 = 



rA[e ] 



2ttvqt — B 



(25) 



Calculating the ladder contribution to the conductivity 
with the renormalized vertex V 2 we get \V2l . 

Finally, we calculate in the ^-approximation the quan- 
tities A[en], B[eo] and C[eo] introduced in ||HJ, ijT^jl and 
ifbfl) . The ^-approximation consists in replacing J2 P ~ * 
v all prefactors before the Green's functions being 
taken at the Fermi level. The pairs of Green's functions 
for the same chirality are the easiest: 



^ ] G u +G u — J 



(£„ - i(u + fi))(£ u - iu>) 



2itv u tI 

flT + 1' 



(26) 

and analogously for the d chirality. The cross-products 
on the other hand depend on the gap A at the Fermi 
surface = ( T A: 



G u +Gd 



2-kvqtI 



Qt + I -z2At' 



(27) 



and the expression with A «-> —A for the crossproduct 
Gd+G u . Hence we find, taking f2r — * and in the band 
of frequencies where 1=1: 



A[e ] 
2tti/qt — B[eo] 



= -2irv Q e' (p F ) 



1 + 4A 2 r 2 ' 



= 2-kvqt 



2A 2 



1 + 4A 2 r 2 ' 



(28) 
(29) 



and the expression lfl6ll for C[eo]. The product 
2ttis 6 / (p f ) = p F , so that aA[e ] = -A 2 r 2 /(1 + 4A 2 t 2 ). 

In summary, we have shown that the vanishing of the 
2D ISHE within the Rashba model of spin-orbit cou- 
pling is specific to both the linear-in-momentum Rashba 
coupling [ill and to the parabolic energy band disper- 
sion. For a general band dispersion, the ISHE within 
the Rashba model is finite, albeit strongly dependent not 
only on the Rashba coupling a but also on the details 
of the band dispersion. In general, the impurity vertex 
corrections suppress the intrinsic cr sH from its universal 



value of e/8ir, reducing it to some small nonuniversal 
finite value of order a 2 . Our result, combined with re- 
cent theoretical demonstrations of finite ISHE in the 3D 
Dresselhaus model [H} and the Luttinger model [2J, es- 
tablish that the intrinsic spin Hall effect is theoretically 
well-defined, but non-universal. 
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